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Abstract 

We have examined the efficiency of baryon-number transport mechanism across the 
phase boundary in a cosmological quark-hadron phase transition through the proper 
estimate of baryon-number penetrability (Eh)- F° r this purpose we have derived first 
the double-pair creation probability Pj, in terms of single-pair creation probabilty per 
unit time and unit volume K m and then obtained an analytical expression for Our 
calculation is free from the uncertainty of the value of double-pair creation probability 
per unit time and unit volume which was used as a free parameter in earlier calcula- 
tions. Finally the variations of double-pair creation probability P^ as well as (S^) with 
temperature are shown and compared with other known results. 
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1 Introduction 



Recently there is a growing interest in the primordial nucleosynthesis in an Universe 
having an inhomogeneous distribution of baryon number. The most popular mechanism 
for the generation of baryon number inhomogeneities has been suggested as a first-order 
QCD phase transition [1-8]. The present Universe is presumed to have undergone sev- 
eral successive phase transitions associated with symmetry breaking in its early stage. 
A QCD phase transition might have occurred at approximately 10~ 6 sec after the big 
bang in the early Universe. Even though many aspects of phase transition are not a 
settled issue, yet, the QCD phase transition is expected to be a first-order one from 
lattice QCD results [9]. Therefore, a large local fluctuation in the baryon to photon 
ratio %/n 7 could arise which might subsequently modify the standard picture of the 
primordial nucleosynthesis (PNS). Such a fluctuation might also result into the forma- 
tion of quark nuggets and/or quark stars etc. [1-8]. The basis for the production of 
isothermal baryon number fluctuation lies in the separation of cosmic phases as follows 
[1,5]. Initially the Universe is in the quark-gluon plasma (QGP) phase at a high temper- 
ature and the net baryon number resides entirely in the quark-gluon plasma phase and 
is distributed homogeneously. As the Universe expands, the temperature drops to the 
critical temperature T c where the quark-gluon plasma exists in thermal and chemical 
equilibrium with the dense and hot hadron gas. Subsequently, the expansion requires 
a continuous conversion of QGP into the hadron phase. The phase transition is com- 
pleted when all the quark-gluon plasma has been converted to hadron phase and all the 
baryon number residing finally in the hadron phase is distributed homogeneously. The 
magnitude of baryon number fluctuation is estimated by the baryon contrast ratio R eq 
of the net equilibrium baryon number density in the quark-gluon plasma (QGP) phase 
to that in the hadron gas (HG) phase, i.e., R eq = nQ GP /n^ G which is evaluated at the 
critical temperature and chemical potential. The baryon number density inhomogeneity 
arising due to such a quark-hadron phase transition will thus alter the yields of the PNS 
which is a function of R eq and the neutron to proton (n/p) ratio [10]. Several theoretical 
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attempts have been recently made to determine the values of R eq [11-14]. 

All aforementioned calculations of R eq assume that the Universe is in complete ther- 
mal and chemical equilibrium during the phase transition [15] . But in reality, the cosmic 
first-order phase transition will necessarily result in deviations from thermal and chemi- 
cal equilibrium because the low temperature hadron phase is not nucleated immediately 
at the critical temperature T c . A generic feature of quantum or thermal nucleation 
theory is that the nucleation rate does not become large unless the temperature has 
dropped below T c [16]. The magnitude of these deviations will depend on the efficiency 
of heat and baryon transport across the transition front during the phase transition. 
Latent heat transport is carried out by the motion of the boundary wall which converts 
volume of one vacuum into another. This vacuum energy difference acts as a source for 
particle creation. However, the latent heat or entropy could also be carried out across 
the phase boundary by neutrinos apart from surface evaporation of hadrons (mostly 
pions). Baryon number is not thermally created in the hadron phase. Therefore, it 
actually flows across the boundary by the strong interaction physics if the hadron phase 
is going to have any net baryon number at all. 

There are two limiting situations governing the efficiency of baryon number transport 
across the phase boundary [5.17]. In the first situation, when the bubble-like regions of 
quark-gluon plasma are small in size, the phase boundaries move slowly compared to 
the baryon number transport rate. Hence the baryon number is quickly and efficiently 
transported across the boundaries, so as to maintain the chemical equilibrium. In the 
second situation, the boundaries of the bubble-like regions move rapidly compared to 
the baryon transport rate so that the net baryon transport process becomes inefficient. 
Therefore, chemical equilibrium may not be achieved on the time scale of the phase coex- 
istence evolution and, consequently, the baryons could be concentrated in the shrinking 
bubble-like regions of the QGP. The magnitude of the resulting baryon number fluctu- 
ation at the time of phase decoupling will still be given by the ratio R' (say) of baryon 
number densities in the two phases but now each phase will have a different baryon 
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chemical potentials. Thus, the final ratio R' will be larger than the equilibrium ratio 

Req- 

It is clear that in both the abovementioned situations isothermal baryon number 
fluctuations will result. However, the magnitude of R eq is in the limit of efficient baryon 
number transport mechanism which is the aforesaid first situation. The efficiency of 
baryon number transport across the phase interface will be given by the bulk properties 
of the phases and the estimate of baryon number penetrability [5,17]. The baryon number 
penetrability from the hadronic side is defined as the probability that a baryon which 
approaches from the hadronic phase will dissociate into quarks and pass over into the 
QGP phase. Similarly from another side, we could define the probability S g that the 
quarks which approaches towards the phase boundary will form a color singlet baryon 
and pass over to the hadron phase. The thermal average of these probabilities are related 
by the detailed balance 

k /^<E ? > = f b _- b (E h ) (1) 

where f q -g (f b _i) is the excess quark flux (baryon flux) in the QGP (HG) phase, respec- 
tively, and the dimensionless quantity k takes values from 1/3 to 1. The value k = 1/3 
implies that baryon number predominantly passes over the phase boundary as preformed 
three-quark clusters in the QGP phase. Similarly k = 1 signifies that baryons will be 
predominantly formed by the leading quarks which cross over the phase boundary into 
the hadron phase. A low value for (Eh) (or (< 1) implies an early drop out of 

chemical equilibrium during the separation of phases and thus might lead to the large 
amplitude of baryon-number fluctuations. 

The estimation of baryon number penetrability has been attempted recently [5,17,18] 
by many authors within the framework of the chromoelectric flux tube (CFT) model. 
Here an energetic leading quark in the QGP is assumed to pick up an antiquark forming 
thereby an expanding CFT whose decay through a single or double pair creation results 
in the formation of mesons or baryons, respectively. Naturally, the calculation of (S/j) 
within CFT model requires the values of the single as well as coherent double pair cre- 
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ation probabilities per unit space per unit time, called K m and Kb, respectively. Although 
K b is a crucial factor for the formation of a baryon yet no theoretical expression for it 
based on strong interaction physics exists in the literature. 

In this context the following points pertaining to a recent paper by Jedamzik and 
Fuller [17, referred to as JF hereafter] become particularly relevant : (i) JF have ex- 
tracted the numerical value of Kb by analyzing empirically the ratio of baryons to mesons 
produced in the e + e _ — > hadrons experiment conducted in the accelerator laboratories, 
(ii) JF allow for the full angular range < 9 < n where 9 is the polar angle of the 
leading quark's velocity with respect to the normal to the phase boundary, (iii) The 
final expression for (E^) written as an integral over the quark distribution functions is 
computed by JF entirely numerically, (iv) Finally, the threshold energy in JF's work 
increases with the temperature in an unreasonable manner. 

The aim of the present paper is to extend/modify the theory of JF [17] in multifold 
respects : (i) By regarding the unconnected double-pair creation event as a succession of 
two single-pair events, we show in Sec. 2 below, that the double-pair creation probability 
Pb (over finite time duration) can be obtained in terms of the K m parameter itself, (ii) 
Since the leading quark's velocity should be directed towards the phase boundary, we 
point out in Sec. 3 that the corresponding polar angle should be restricted to the range 
< 9 < 7r/2. (iii) We derive in Sec. 3 an analytical expression for the thermal average of 
the baryon number penetrability so that the dependence of (Eh) on the relevant variables 
becomes more transparent, (iv) We show in Sec. 4 the comparison between our and JF 
models by plotting P b as well as (E h ) as functions of the temperature when the quarks 
are assigned either the constituent or current mass, (v) Finally, we carefully analyze the 
temperature-dependence of the threshold energy in Sec. 4 to obtain an elegant formula 
for (E h ) in the high T limit. 
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2 Decay Statistics of Flux Tubes 



Chromoelectric flux tube (CFT) models provide us a phenomenological description re- 
garding the formation of a hadron from quarks and gluons [19]. These models assume 
the existence of a chromoelectric field between two oppositely coloured quarks where the 
chromoelectric field strength is assumed to be constant in magnitude and independent 
of the separation of the quarks. These fields can be thought of as being confined to a 
tube of constant width known as flux tube. Lattice QCD results justify the existence of 
such a flux tube [21]. Initially, chromoelectric flux tube models were used to explain the 
processes such as e + e~ — > hadrons [22,23]. Later these models have also been used to 
describe the spectrum of mesonic and baryonic resonances [22,24]. Recently flux tube 
models have been employed to estimate the meson evaporation rates from quark-gluon 
plasma produced in relativistic heavy-ion collisions [20] and also to estimate baryon- 
number penetrability across the phase boundary in cosmological quark-hadron phase 
transition [17,18]. 

If the CFT is regarded as an unstable system then its stochastic properties can be 
conveniently discussed by introducing the following abbreviations : SP = survival prob- 
ability, DP = decay probability, m = mesonic or single-pair production channel, b = 
baryonic or connected double-pair production channel, b* = disconnected double-pair 
production channel. Let us now take up the the earlier model used by Jedamzik and 
Fuller [17] and a new proposal due to us. 



2.1 Jedamzik and Fuller (JF) Model 

In analogy with QED, at a perturbative field theoretic level, the JF model amounts to 
connected diagrams shown in Figs. l(m) and 1(b). The corresponding DP's per unit 
time per unit volume are denoted by K m and Kb, respectively, having a total k = K m + k^. 
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Figure 1: (m), (b), (b*) : Diagrams drawn in analogy with perturbative QED, repre- 
senting decay into single, connected double, and unconnected double qq pair channels, 
respectively. The crosses stand for the chromoelectric field as a source. The dotted lines 
denote gluons. 

The space-time integrated parameters are called as 

rto 

W m = Km V(t)dt (2) 

JO 



and w = w m + Wb, V(t) is the instantaneous volume of the flux tube and t ~ Eja is 
the maximum time upto which it expands when the incident quark has energy E. From 
the conservation of energy and momentum parallel to the phase boundary, the above 
integral could be written as in [17] 



with 6 being the angle of incidence of the leading quark relative to the normal to the 
phase boundary, A, the radius of the flux tube and a the string tension of the flux tube. 




(3) 




(4) 



Then, according to the classical theroy of radioactivity, the net SP and DP counting all 
channels would be 

P s = e~ w (5) 
Pa = l-e- w (6) 

Since the relative weights for the m and b channels are w m /w and Wb/w, respectively, 
the corresponding channelwise DP's are as used by JF : 
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1 — exp 

Although be estimated from Schwinger's single-pair production mechanism [25] 

applied to QCD, nothing a priori is known about k^. However, JF have extracted 
empirically the magnitude of K h from the analysis of e + e~ — > hadrons data. 

2.2 Our Proposal 

Suppose we retain information about the single qq pair production of Fig. l(m) but 
ignore the diagram 1(b). Then, the double-pair production event may be looked upon as 
a sequence of two unconnected single-pair creations within time to as shown in Fig. 1(6*). 
Since the flux tube is now characterized only by the parameter w m , hence radioactivity 
theory would give the following expressions for the net SP and DP over the duration 
time to : 

P* = e~ Wm (9) 

oo 

p; = i-e-^ = £K (10) 

n=l 

where the DP of having exactly n successive single-pair events is a Poissonian viz. 



Therefore, in our approach the DP for the single and double- pair would be obtained 
from 

K=Pl = e~ Wm w m (12) 

P* = p * 2 = e ~ w - (13) 

In contrast to the JF approach (7,8), our proposal involves only one parameter w m (or 
K m ) as an input. 



3 Analytical Estimation of Baryon number Penetra- 
bility 

The main result of JF's work is contained in their baryon penetrability equation which 
gives the thermal average of the baryon penetrability integrated over the qq distribution 
function. We wish to make two plausible comments on this result. Firstly, the angular 
integration range in their calculation < 9 < n may not be appropriate because if the 
leading quark is to enter the phase boundary then it should have < 9 < ir/2. Therefore 
we suggest that the corrected expression should be 

W = f ^ d9 r dE ±1(9) P b (E, 9) (14) 

Jb-b Jo jE th dEd9 

Here f b _i, the excess baryon flux in the hadron phase, is given by [17] 

fb-b = % £- 2 K + T) T 2 exp(-m b /T) (15) 

where mj, jj, b and g b are the baryon mass, baryon chemical potential and degeneracy fac- 
tors, respectively and x\ is the component of the leading quark velocity perpendicular 
to the phase boundary, and E t h is the threshold energy. 

Secondly, JF have done their subsequent calculations numerically so that the dependence 
of (E/j) on the relevant parameters remains somewhat obscure. We suggest that an ana- 
lytical estimate for the seemingly complicated expression (14) is very desirable because 
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that would make the dependence of (E^) on the relevant parameters more transparent. 
For this purpose we proceed as follows : 

Using Boltzmann approximation(E/T 1), the differential excess quark number density 
in a given energy interval dE and in a given interval of incident angles d6 becomes 

dn n 



"q-q 



% E 2 exp(-E/T) sine (16) 



dEdO T 2tt 2 

where /i q , g q are the quark chemical potential, and the statistical weight of quarks, 

respectively. 

Then Eq.(14) becomes 

(E/,)«C ds s dE E 2 exp(—E/T) P b (E,9) (17) 

J0 JE th 

where 

c • = a » t (18) 

Due to the rapidly damped factor exp(—E/T) in Eq.(17) most contribution to the en- 
ergy integration comes from around the threshold energy E t ^. Therefore, it is convenient 
to make the transformation 

P = ; E = E th {l + -^p} (19) 

Upon using the identity / °° dp exp(-p) = 1, Eq.(17) reduces to 

(£/,) » C f'dssT E 2 h exp(-^) P b (E th , 9) (20) 
jo 1 

within correction terms of order T/E th . Now, at large angles 9 — > 7r/2, s — > 0, E th — > oo 
implying that exp(—E th /T) gets heavily damped. Hence, most contribution to the 
angular integral must come from around the forward direction 9 — > 0, s — > 1. For 
convenience define 

E t h = E th \ e=0 

= m. + Bn- 1 (21) 
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A 



: P h (E th ,6 = 0) 

E t h — E t h 
T 



(22) 



(23) 



Here B is the bag constant and n q is the quark plus antiquark density. The above 
eq.(21) is the threshold energy condition in the forward angle direction (6 = 0) and the 
next eq.(22) is the corresponding probability. Using the identity f™ b ^ T d\ exp(-X) m 1, 
eq.(20) yields 



" 2 



C T " Eth exp(- E th /T) A {1 + 0{— )} 



Substituting Eq.(18) for C, we arrive at the desired analytical estimate 



\H 9q I E th \ 

lib 9b \m b J 
The probability function in JF model 



T 



e -{E th -m b )/T p & | 1 + (_)| 



P h = 



exp{-b E th } 



can be replaced in our model as 



with 



a = 



p: = 



K b 
K m + K b 



Wt 



b = ( Km + Kb ) 



ttA 2 
2a 2 



W r . 



K r . 



2a 2 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



Here the string tension of the flux tube a ~ 0.177GeV 2 , and single-pair creation proba- 
bility per unit time per unit volume K m is given by [25] 



Km. 



47T 3 



exp(— Tim 2 /a 



(30) 



where m q is the quark mass. Since there is an ambiguity in choosing the precise value 
of m q , we allow in the sequel both the possibilities viz. m q = 300MeV (constituent 
quark mass) and m q = lOMeV (current quark mass). Eqs.(25), (26), (27) form the 
main algebraic results of the present paper. 
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4 Results and Discussions 



In Fig. 2, we have shown the variation of double pair creation probability at threshold 
in the forward direction with the temperature for a small value of b = 0.047 GeV~ 2 
corresponding to the constituent quark mass (m q = 330MeV). Here we have taken the 
value of the ratio a(« K m / K b ) rs 1/5 as in [17] 
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Figure 2: Variation of double pair creation probability P b (P b *) at threshold in the 
forward direction with the temperature T when the quark has constituent mass. The 
dashed curve represents our model calculation (cf. Eqs.(13), (27)) whereas the solid line 
represents the JF model (cf. Eqs.(8), (26)). For remaining symbols see Eqs. (28,29,30) 

Fig. 3 is also the same as Fig.2 except that a higher value of b = 0.32 GeV~ 2 is taken 
corresponding to the current quark mass (m q = lOMeV) and the ratio a ^ 1/20. 

Clearly, the difference between the predicted values of P b in our model and JF model 
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Figure 3: Same as Fig.2 except that current quark mass has been chosen in both 
models. 

must be attributed to their different premises. Our approach has the advantage of 
working with only one decay parameter K m but ignores the connected diagram 1(b). 
The JF approach has the advantage of including the graph 1(b) but at the cost of 
bringing-in the additional unknown decay parameter Kb- 

Next, we turn to the calculation of the baryon number penetrability (T, h ) based on 
the analytical result (25). The probability P b can have either the JF form (26) or our 
form (27) with the remaining parameters having been set as for u, d flavours in QGP 
sector and nucleons in the hadronic sector. Figs. 4 and 5 display the dependence of (E^) 
on the temperature when the quark has the constituent and current masses, respectively. 

In Fig. 4, when the quark has constituent mass, our curve lies below that of the JF 
model. Our small predicted value (E h ) ~ 0(1CT 3 ), incidentally, justifies the results 
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Figure 4: Variation of baryon number penetrability (E^) (cf. Eqs.(25), (26), (27)) with 
the temperature T when the constituent quark mass is chosen. Other notations are the 
same as in Fig. 2 

obtained by Sumiyoshi et al.[18] and Fuller et al.[5]. However, in the opposite limit, 
i.e., for current quark mass, Fig.5 reveals that our graph is consistently above that of 
JF. Our predicted higher value of (E^) ~ O(10 -2 ) seems to agree with the original 
result obtained by Jedamzik and Fuller [17] using numerical quadrature involving the 
unphysical angular range < 9 < ir. 

In passing we note that the typical time upto which the CFT expands is t ~ E th /a ~ 
10~ 24 s which seems to be one order of magnitude smaller than typical nuclear time scale 
of 10~ 23 s. Can the b* channel of Fig.l be important under such a circumstance? To 
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Figure 5: Same as fig. 4 except that the current quark mass is chosen, 
answer this question we borrow the w m values from Figs. 2-5 and look at the ratio 

H = ^ ^ 2% to 15% (31) 

p* o v ' 

which is sizable. Therefore, our view of regarding the double-pair production event as a 
succession of two single-pair events can be justified. 

Towards the end we wish to comment on a couple of factors on which baryon number 
penetrability crucially depends within the chromoelectric flux tube approach. One factor 
is the choice of the quark mass as seen from Figs. 4 and 5. Another factor, in which 
(Eft) is very sensitive, is the threshold energy E t h in the forward direction. Following 
the idea of Jedamzik and Fuller [17], this threshold energy E th consists of two parts : 
one is simply the rest mass of the baryon and the other is the interaction energy 
Bn~ l for each quark and antiquark which resides in the QGP. JF have parametrized the 
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interaction energy as ~ 3.7T which apparently grows with the temperature. However, 
in our opinion, the function Bn^ 1 should decrease with T like T~ 3 for a fixed value of 
the bag constant B since the total quark density n q is an increasing function of T. So 
E th would tend to m b in Eq.(25) which yields the elegant high-temperature behaviour 

(£„) ^ A (32) 

remembering that g q = 12, p b = 4 and /x^ = 3/i g . 
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